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Abstract 

J> ' We develop an approach for calculating matrix elements of meson exchange 

currents between 3N basis states in (j J)-coupling and a 3N bound state. The con- 
tribution generated by ir— and p-exchange are included in the consideration. The 
matrix elements are expressed in terms of multiple integrals in the momentum space. 
We apply a technique of the partial wave decompositions and carry out some angu- 
. lar integrations in closed form. Different representations appropriate for numerical 

calculations are derived for the matrix elements of interest. The momentum depen- 
dences of the matrix elements are studied and benchmark results are presented. The 
approach developed is of interest for the investigations of deuteron-proton radiative 
^5 | capture and 3 He photo- and electrodisintegration when the interaction in the initial 

or final nuclear states is taken into account by solving the Faddeev equations. 



^ ! 1. Introduction 

H , 

Various approaches have been developed to calculate the amplitudes of photo- and electro- 
disintegration of 3 H and 3 He and the radiative capture of nucleons by deuterons. Methods 
employed and results obtained recently in studying these reactions can be found in |l||-[|T8||. 



Review article [HI provides additional important references. 



A common feature of investigations [l|]-[|T7] is the use of the solutions of the Faddeev 



type equations for the 3N bound states. While the results of U have been obtained 
with the s-wave Malfliet-Tjon NN interaction, realistic models of nuclear forces (e.g., the 
Argonne, Bonn, Paris and Reid soft core potentials) have been employed in p|-||17||. 



In H, 0, f§, and JTT|-JT8[ the rescattering effects (the interactions in the initial or 
final states) are taken into account in different ways. For example, in |]TT| the wavefunction 
(WF) of the pd-system for the reaction 

d + p^ 3 He + 7 (1.1) 
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has been calculated partially in the coordinate and partially in the momentum space 
representation applying the method of continued fractions, whereas in JL6] the initial state 
interaction for the same reaction has been considered applying an iterative procedure to 
a Faddeev-like equation in the momentum space and summing up the resulting multiple 
scattering series via the Pade method. 

Calculations and |17| are based on separable approximation [ll|,[20| constructed 



with the help of the Ernst-Shakin-Thaler method^TJ. The use of correlated hyperspher- 
ical harmonic WFs for the 3N systems is mentioned in 



An essential element of all these treatments is a model of the electromagnetic (EM) 
current for the 3N system. In JjJ, P|- [[PT|| the current is chosen following the basic idea of 
the impulse approximation, i.e., it is assumed that the current is an one-body operator. 
Two-body contributions to the current (interaction or meson exchange currents (MEC)) 
have been included either implicitly [|l|], |8[|-|tL7|| as pointed out below or in an explicit 
form [IS]. Earlier the contributions of MEC to the magnetic multipoles have been 

considered in investigations [ 22 ] — [ 23 1 of the radiative capture ( |1 . 1|) . 

In this connection, it is pertinent to recall (see, e.g., |E3],[E3]) that the necessity of the 



;i.2) 



inclusion of MEC results from the continuity equation (CE) 



o 



or its equivalent 



H,f(Q)} = \P,j(0) 



;i.3) 



where j' 



j(x)) is the current density operator and P is the total momentum 



operator for a given system (nucleus). The nuclear Hamiltonian H = K + V consists of 
the kinetic energy part K and the interaction part V . The latter is isospin dependent, 
and, in general, can be nonlocal, i.e., dependent not only on the nucleon coordinates but 
also on their derivatives. 

For a system with such a Hamiltonian where the forces can depend arbitrarily on 
nucleon velocities a conserved EM current can be constructed following, for example, a 



prescription [26 



As is well known (see, e.g., ||27||,|pE|), the CE for the current operator is not a sufficient 
condition for ensuring gauge independence (GI) of the reaction amplitudes. Indeed, the 
GI requirement reads 

<f (/|i/i(0) \i) = 0, (1.4) 

where g M is the four momentum transfer. Eq.( |1.4j ) is fulfilled if the operator j M (0) meets 
Eq. (|1.3|) and both the initial \i) and final [/) nuclear states are exact eigenvectors of 
H. The implementation of these requirements is a quite demanding task in practical 
calculations because one has to take care of the definite consistency between the current 
and the Hamiltonian and "good" WFs are needed both for initial and final states. In this 
respect, an exception is the calculation of the charge form factor (FF) for a nucleus (for 
instance 3 He ). In this case one can show that the one-body (j^(x)) and many-body 



x)) contributions to the decomposition j(x) = j^(x) + jl meson ](x) satisfy 
see Refs.l 



[31]). Of course, it does not mean that one can avoid 



meson (j\- mei 
Eq. (|1.4|) separately 
the inclusion of MEC. 

For the nuclear processes with real or virtual photons Eq.([OJ) is useful for an effective 
account for MEC effects. Using Eq.( |1.4| ) the matrix elements of the longitudinal compo- 
nent of the EM current can be expressed through the corresponding matrix element of its 
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time component. This replacement makes sense due to the fact that the latter is subject 
to the meson exchange influence to a lesser extent than the space component of the cur- 
rent (cf. classification of the MEC given by Friar |32|] and discussion in Sect. 2). This 
trick is often employed in the treatments of inelastic electron scattering on nuclei in order 
to express the differential cross sections in terms of the respective structure functions 

1,0-0. 



12 -15 , 33 



One can also write the amplitudes in a GI form via extended Siegert theorem ]34 



36] expressing them through the strengths of the electric and magnetic fields and the 



generalized electric and magnetic dipole moments of a system (nucleus) 

The Siegert theorem is widely used as a recipe to incorporate a part of MEC effects 
within the long wavelength approximation. However, all the applications of the Siegert 
theorem imply that the initial and final states involved in the transition matrix elements 
are solutions of the Schrodinger equation with the same Hamiltonian. This important 
point has been violated in some previous calculations (see, e.g., |10j and [p7|l ). 

The initial and final nuclear states have been treated consistently in studies 0, @, @- 

rm 



and 



of inelastic electron scattering on 3N nuclei and the radiative capture. 
The methods developed in three-particle scattering theory |p8|-|i0| provide a rigorous 
description of the nuclear states belonging to the continuum. All relevant components 
of realistic NN potentials can be included. In |12|— 116| rescatterings in the two- and 
three-body channels for the 3 He breakup by electrons are treated on the same footing by 
solving the Faddeev-like integral equation 



\U K ) = tG (l + P)j K {Q) Abound) + tG P \U K ) 



1.5) 



where j K {Q) {k = 0, ±1) are the spherical components of the Fourier transform of the 
current operator. The calculations have been performed in the pga-representation, where 
p and q are the Jacobi momenta and a stands for discrete quantum numbers in (jJ)- or 
(LS)-coupling. As is noted in Hl4| , the block j^{Q) Inbound) is the only essentially new 



structure in comparison to the nucleon-deuteron scattering theory. 

Only the one-nucleon part of the EM current has been utilized in ||I2"|-||16||. Unlike 
this, in studies |§-[0], |41| of the two-body 3 He disintegration by real and virtual photons 
at the intermediate energies special attention has been paid to the MEC contributions 
with no final state interaction included. The two-body currents ||2], [J3| due to the n- and 
p-meson exchanges and MEC with the excitation of the A-isobar in intermediate states 
have been taken into account. In 0-0, pl[ the reaction amplitudes have been evaluated 
with a convenient parameterization [44j] for the Faddeev components of the 3 He WF. 



This component has been computed with the Reid soft core potentials in [|45 |. 

The simultaneous consideration of the MEC and rescattering effects in the theory of 
the reactions still remains to be carried out. The purpose of the present paper is to 
develop a method for the calculation of the matrix elements of the two body current 
between the 3N basis states and the 3N bound state. A particular emphasis is laid on a 
compatibility with the approach implemented in |T^[-|[TB[. We treat the matrix elements 
in the momentum representation and use the basis in the space of 3N states formed by 
the vectors with the quantum numbers in (j J)-coupling. 

We aim to elaborate a technique which can allow us to provide the GI of the reaction 
amplitudes in the calculations with the WFs for realistic models of nuclear forces. 

This paper is organized as follows. In Section 2 the transition matrix elements of the 
two-body currents are transformed to forms which are suitable for numerical calculations. 



3 



They are expressed in terms of two- and four-dimensional angular integrals. Section 3 and 
Appendices A-C are devoted to the reduction of the angular integrals. The partial wave 
decomposition of the p-meson exchange currents is given in Appendix D. We show two 
ways how one can analytically carry out three integrations over the angles. The explicit 
expressions for the matrix elements of interest are presented. Section 4 gives benchmark 
results in the case of the 7r-meson exchange currents. We use numerical methods to eval- 
uate angular integrals employing MEC matrix between the 3N basis states and compare 
the results with ones obtained within the framework of the approaches proposed above. 
Section 5 contains the summary and a discussion. The Appendices contain detailed ex- 
pressions for the matrix elements and the integrals involved. In particular, in Appendix 
D some angular integrals appearing in the derivation are calculated in a closed form. 



2. Contributions of Two-Body Currents 
Amplitudes of the reactions 



can be written as 



7 + 3 He -> p + d, e + 3 He -> e' + p + d, 

7 + 3 He — > p + p + n, e + 3 He — > e' + p + p + n 



A''" = 



f(Q)%), = 0,1,2,3), (2.r 



where j^{Q) is the Fourier transform of the current density ), the initial and 
final ^ states are eigenvectors of one and the same Hamiltonian H. 

Our consideration refers to the case where the current j(Q) = j^(Q) + j' 2 '(Q); i- e -; 
consists of a one-body current 



3 [1 \Q) = j(Q; i) + 3(Q\ 2 ) + KQ; 3) (2.2) 

and a two-body current 

3 [2 \Q) = j(Q] 1, 2) + j(Q; 1, 3) + j(Q; 2, 3). (2.3) 

The operators j^(Q;a) and j^(Q;b,c) depend on the variables of the nucleons with 
the labels a and b, c, respectively. Taking into account these decompositions and the 
permutation properties of the states, one can write 



f(Q;l)+j»(Q;2,3) % ) . (2.4) 



j\P = 3 (y { f ] 

The method for calculation of the one-body contributions to amplitude ( |2.4|) has been 
worked out earlier and used in studying the elastic electron scattering off 3 He |29|| , the 
3 He disintegration by electrons [0-[|15| and the proton-deuteron radiative capture [|16|]. 
The detailed description of the corresponding formalism can be found in the above articles. 
Here we concentrate on a treatment of the two-body contributions. 

In the contemporary theories of EM interactions with nuclei, where many-body cur- 
rents are generated by meson exchange, the time component jg 2 ' is of order v 2 /c 2 com- 
pared to j'o^ which is of order (v/c)°, while j ' 2 ' ~ j M = 0(v/c). Therefore, neglecting 
the v 2 /c 2 -contributions, we assume 

Jo - jf 1 (2.5) 
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and 



J =J [1] +J [2] - 



(2.6) 

In other words, the dominant MEC contribution to the reaction amplitude originates from 



the space part of the current, ' j(Q; 2, 3) 

To get the matrix elements of interest let us project the vector j^(Q; 2, 3) |^) on the 
basis 

\pqa) = \pq (ls)j (A|)I (jl)JM; [t\)TM T ) i (2.7) 

used throughout Refs.||l|-[||], ||, @-||. 

In Eq. ( |2.7| ) the subscript "1" is to indicate the definite choice of the Jacobi momenta 



P= -(^2 - k 3 ), 



q = -(2/ci - k 2 - h) 



{21 



and the angular momentum coupling scheme as it is depicted in Fig.|l]. The nucleon mo- 
menta are denoted by k a (a = 1, 2, 3). The orbital angular momentum, the total spin, and 



3Q 



P / 



©1 



*©2 



Figure 1: Jacobi momenta. 



the angular momentum in the two-body subsystem 2-3 are /, s and j, respectively. The 
orbital angular momentum A of nucleon 1 and its spin couple to the angular momentum 
/ of this particle. The angular momentum of the 3N system and its magnetic quantum 
number are denoted by J and M. The total isospin of particles 2 and 3 is t, while T, Mt 
refer to the isospin of the 3N state and its projection. The set of states ( |2.7|) is orthonor- 

*$r } ) and 
by requirement 



malized and complete (see, e.g. fl40| , Section III). Since the nuclear states 
are antisymmetric we can consider the subset of vectors (|2.7|) constrained 



/ + s + t = 2n + l, 



(n = 0,l,2,...). 



It is convenient to split the whole array of discrete quantum numbers a, separating 
the isospin part of vectors 



where 



\pqa) = \pqotj) \a T ) 



\aj)= (ls)j,(X^)I,{jI)JM 



(2.9) 
(2.10) 
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and 



\a T ) 



(t-)TM T 



(2.11) 



The matrix elements of the two-body currents we shall consider have the form 



k 2 k 3 



J(Q; 2,3)| k 2 k 3 ) = 6(k' 2 + k' 3 -Q-k 2 ~ k 3 ) j (k' 2 - k 2 , k' 3 - k 3l 2,3) , 



f2.12) 



or in 3-particle space (without the overall momentum conserving delta-function) 



p'qi 



j[Q;2,3 



pq 



s (<f - q + <5/3 

xj{Q/2 + p'-p,Q/2-p' + p; 2,3 
= S (p 2 + p 3 - <5) j(P2,P 3 ', 2, 3) 



(2-13) 



with 



p 2 = k' 2 - k 2 = Q/2 +p' - p 
p 3 = k' 3 - k 3 = Q/2 -p' + p 

The corresponding matrix elements of the one-body current are 

[k'x \M 1)| h) = 5(k[ - g - Q)M, Q; 1) 
or in 3-particle space and in the lab frame 



p'q' 



M 1) | PQ) = S(p~ p')5{q> - q- ^Q)j(q> + Q\ 1) 



(2-14) 



(2.15) 



(2.16) 



Eqs.( p.l^ ) and (|2.16|) exhibit the important shift prescriptions 5 \ q' — q + Q/oj and S(q'— 

q — |Q). The one-body current matrix element has been dealt with in [p9| , |14ll in great 
detail and therefore we shall concentrate now only on ( p,13|) . 

Because the factor j [Q/2 + p ' — p, Q/2 — p ' + p; 2, 3) in ( |2.13|) has no dependence 
on q we can treat the p— and q— spaces separately. We regard the partial wave decom- 
position of the two-body current formed in terms of the basis states (|2.7| ) or ( |2.9|) . 



p'q' a 



J(Q;2,3) 



pqa 



d£ I di j dpi J dqi ~£ CU'I'f, pi, M' - ,l')C(j, I, J, V, M - t i.) 

Pi q\ 6 



(2.17) 



j(Q/2+p-['-p-[, Q/2-p^ + p^ 2,3 

/ m *(pi -p) s (qi - q) 

xyj^{pi>yi,M-^qv- 



OtT 



Pi' 



91' 



The momentum states \p\qi) and \p'iq\) are formed out of Jacobi momenta of the type 
~M. We obtain 



p'q'a' 



J(Q;2,3) 



pqa 



= C(j'I'J', n', M' - n')C(j, I, J, /i, M - //) 
J dp[ J dpxy*,^, (p[) (a T > \j (Q/2 + p'p[ -pp u Q/2- p'p\ + ppi, 2, 3) | a T ) y j/x (pi) 

• n %-|g'9i + |Ql) 



1 



q* ' • - - 3 

E C(j'l'J', //, M' - fi')C(j, I, J,n,M- li) 



with 



xl 2 {p',p, Q; (l's')j'/j'a T ,, (ls)j/ja T ) 
xl 3 (g', g, Q; (A^)I'M' - (A^)/M - 

Q; {l's')fn'a T >, (ls)jfi, a T ) = 
J dp' J dpy^,{p')(a T , \J(Q/2+p'-p, Q/2-p' + p; 2, 3) \ a T ) 

W, g, Q; (A^)I'M' - //, (A^)IM - //) 

Finally 3^>(p) and 3^,i/(<?) are defined as 

y^ip) = ^C{lsj,ii-m,m,ii)Yi^ m {p)\sm) 



(2.18) 



and 



(2.19) 



and 



yiu{q) = *£,C{\-I,v - m,m,u)Yx,u-m(q)\^m). 



Now we use the result ( |2.19| ) to evaluate the decisive building block (pqa j(Q; 2, 3) 
We obtain 



(2.20) 
(2.21) 

^ bound 



I 1 1 

pqa 



j(Q;2,3) 



= E / P 2 dpq 2 dq (p'q'a' \j(Q; 2, 3)\pqa) (pqa\V bound ) (2.22) 
E / V 2 dp E Cifl'J', li, M' - fi')C(j, I, J,»,M- li) 



xl 2 (p',p, Q; (l's')j'p,'a T ,, (ls)jp>a T ) 
xl 3 (p, q', Q; (A'i) I'M' - //, (\ l -)IM - 



with I 3 given as 



h(p, q\ Q\ (A'~)I'M' - //, (A~)IM -n) = 
dq'y*r, M >-M){^ $ + \QH^ bounty i, M-^q' + ^Q). 



(2.23) 



The angular integration over q' can easily be performed as demonstrated in Appendix A. 
The integral I 2 is more complicated because the operator j also acts on the spin space. 
To proceed we need some specific expressions of the two-body currents. 

Let us return to integrals ( [2.18 ) and examine the spin-isospin structure of the inte- 



grands. In this section we consider model [42[,||43l of 7r-meson exchange currents and 
assume that 

jl [2] ~jseagull _|_ ~jpionic 

The corresponding contributions to j(p 2 ,P3', 2, 3) read 



^2.24) 



j sea9ull (p 2 ,p 3 ; 2,3) 
= F?(v v (p2) (a(2) • p 2 ) <r(3) - v w (p 3 ) a (2) (a(3) • p 3 ) 
x i [7(2) x f(3)] z , 



(2.25) 



■ptonic 



(P2,Ps] 2,3) 



= FY{n-n) (<?(2)-p 2 ) (5(3) -ps) 
x i[f(2)xf(%, 



Vn(P3) ~ Vn(p2 
P2-P3 



where 



v n (k) 



ft 



NN 



2TC 2 m 2 r m 2 + k 2 



tNN 



[2.2Q) 



(2.27) 



further F{ = Ff — F™ is the isovector nucleon form factor, f^NN is the pseudovector ttN 
coupling constant, m n is the pion mass and F n NN{k 2 ) is the strong 7iNN form factor. 
The operator j (p 2 ,P3] 2, 3) can be written as 



j (P2,Ps;2,3) 



F\ EO kK (p 2 ,p 3 ) {<r(2)®a(3)} kK 1 \r{2) x r(3)] 



f 2.281 



where {cx(2) ® °"(3)} fcK is the irreducible tensor product of rank k (k=0,l and 2). 
According to Q2.24j) - (|2.26|) the functions O kK (^2,^3) are 



O 



kK 



^seagull ' ^ •pionici 



where 



O seagull iP2,P?i) 



Y.c{iik^c*){-if v 7r (p 2 )( P2 u+ (-i) fe+ V(p 3 )(P3)- f 



^2.29) 



(2.30) 



with (, £ = 0, ±1 and 



&$LiJfii,Pa) = (P2-P3) (p 2 2 



Kfej) - V n (p 2 )) {p 2 ®P 3 } 



A: re 



P2-P3 

— 9 7 

P2-P3 



K(pa) - ^(p 2 )) E f C(llfc, f , K - f , kO(-) k Go 2 )- € (p3)?-k 



(2.31) 
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Inserting ( |2.28| ) into ( |2.18| ) we get 



h{p,p',Q] {l's')j'n'a T ,, (ls)jna T ) 

= C{l's'j\m\m' B yL')C{lsj,mim a yL) 

m!,m' s mims 

xJ2Mf kK (p',p; Q) (s'm' s \{a(2) ® <r(3)} J sm s ) 



x «|i[r(2) x r(3)y a T > 



The angular integrals 



Mj'^' kK (p',p;Q) 



J dp J dp' Y* m ,(p')d kK (Q/2+pi-p, Q/2-p' + p)Y lm (P) 



are analyzed in the subsequent sections. 

The spin and isospin matrix elements have the form 



(2.32) 



(2.33) 



'm! | {cr(2) ® cr(3)} fcK | sm) = Q\fi¥'C(kss' ', Kinm!) 



2 2 

1 1 k 



{a' T \i[r(2) x r(3)] z \a T ) 
= l2^{-l) 1+t+1 2 +T 'Vtt'C{l\T'^M T M 7 



1 t i! 

1 T' - 

2 2 



1 1 1 

1 i t 

2 2 1 

lit' 
2 2 fc 



(2.34) 



(2.35) 



where a — 2a + 1. It has been used in Eq. ( |2.35 ) that the isospin of the initial state 
T = 1/2. 

Below we develop and study the potentialities of diverse approaches for evaluation of 
the angular integrals ( |2.33| ). Two different method for further reduction of the matrix 
elements are presented in Section 3 and Appendices B and C. While Section 3 and Ap- 
pendix B deal with the 7r-meson exchange currents, the partial wave decomposition of 



the p-meson contributions []42|,||43|| is given in Appendix C. 

Both approaches based on the analytical methods end up with two-dimensional an- 
gular integrations. These integrals and the four-fold ones in ( |2.33|) ( called 4FI in the 
following ) have been calculated numerically. 



3. Reduction of the Multiple Integrals 

In this section we show how one can reduce the dimensionality of the angular integrals 
( [2.33D and transform them to the form appropriate for the numerical evaluation. It is 
convenient to consider the cyclic components 

<r; fe V, p;Q) = JdpJ dp' r; m , (p ') of (p 2 ,p 3 ) Y lm (p) (3.1) 

of tensors ( |2.33|) . We denote 

Pn = lQ + (-l) n (p'-p), (n = 2,3). (3.2) 
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As it follows from Eq.(ggg) 

M^(p',p;Q) = Mt; k «(p',p;Q; seagull) 

+ M^ kK (p',p;Q;pionic). (3.3) 

From Eqs.(CT),(^3TD we get 

M imil kK (p'>P', Q\ seagull) = (-l) fc+K +^(ll£;, u, -«) 
x (|Q„ (G^'(2) + (-l) fe+1 G^'(3)j + G|^(2) + (-l) fc G^(3); 



(3.4) 



and 



where 



M i ;^ kK {p\p-Q- pionic) = (-l) fc+ ^(llfc,^ 3 ,-K 



K2K3 

1/n /n Tpl'm' /n Tpl'm' _i_ /n Tpl'm' nzpl'm' 



(3.5) 



and 



G&'W = JdpJdp'Y i r ml (p')Y lm (p)v( Pn ), (3.6) 

Gfir» = Jdpjdp'Y; ml (p')d,Y lm (p)v(p n ), (3.7) 
(n = 2,3) 

^ = Jdpjdp'Y i r ml (p')d lll Y lm (p)H(p 2 ,p 3 ), (3.8) 

^S'^ = fdpJdp'Y* m ,(p f )d^ 2 Y lm (p)H(p 2 ,p 3 ), (3.9) 

7 /i^ 2 M3 = / ^/ (P d^d^ (p) H(p 2 ,p 3 ). (3.10) 



The explicit expressions for the scalar functions f (p n ) and H(j>2,p3) are given in [Appendix 
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We direct the z-axis along the vector Q. The transformation properites of the inte- 
grands in Eqs.( |3.6| )-( |3.10| ) under rotaions about the z-axis allow one to conclude that 

f^il'm' — X Tpl'm! _ j: 

^Im i"J u m',mi c Imfii u m',m+fii 

^lmfi( n ) = 3m',m+fi, ^/mpti(U2 = ^m',m+A*l+A»2 ) (^-11) 

Tpl'm' r 

- r Zm/xi/^2M3 ~~ y m',m+/ii+/i 2 +M3- 

We express the components of the vector d = p'—p through the spherical harmonics 
with the help of the relation p K = \p | Yi K (# p , P ). Separating the depedence on the 

azimuthal angle we write the spherical harmonics in terms of the normalized associated 
Legendre polynomials Y im (9,<fi) = \ j\p2m exp(zm0) P™(cos^). 

Being invariant under rotation of the coordinate system, the functions v(p n ) and 
H{j>2,pz) depend on p- Q, p ' ■ Q and p-p ' or equivalently on p, p', Q, x, x' and cos(0 — (j)'), 
where x and x' are the cosines of the polar angles of the vectors p and p ' . 

So, we arrive at the integrals of the type (we display here the most complicated one) 

+1 +1 

JdxJ dx'pf(x')pr(x) (312) 

x Pf (x x )Pr (x 2 )Pf (x 3 )R M (p, p', Q; x, x') 
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with 



2tt 



2- 



R M {p,p',Q;x,x') = J d<f) J d(j)' 



x exp[zM(0 — (f)')}f(p,p', Q; x, x', cos(0 — (j)')). 



(3.13) 



In formulae ( |3.12 ) and ( 3.1 3| ) we denote a linear combination of m', m, \L\, fi 2 , ^ by M, 
the variables X\,x 2 and £3 should be identified with x or x' . 

Introducing the new integration variables t = — 0' and t' = + 0' we get 



2tt 



R M (p,p',Q;x,x') =2 I dt(2ir-t)cos(Mt)f(p,p',Q;x,x',cost). 



(3.14) 



The explicit expressions for the functions F and G in the form of Eq. ( |3.12|) can be 
found in [Appendix D. 



It should be noted that the above procedure can be also used in the case of other 
models of MEC. The next step depends on the specific form of the current operator. 
For MEC given by ( |2.24| )-( [2.27| ) and the strong form factor chosen in the form 

(3.15) 



(3.16) 



F-kNN (k 2 ) = (AJ - ml) j (AJ + 



the function f(p,p', Q; x, x', cost) can be decomposed into a sum of pole terms 



\Q±{p-p')] +fi 



(a + boost) n , 



(n = 1 and 2), 



where \x = or A n . In its turn, function (|3.14 ) can be written (see Apendix E for 
details) ClS db superposition of the integrals 



2tt 



V(a, b, m, n) 



dt (2ir — t) (a + b cos t) n cos (mt) 



(3.17) 



calculated for integer n and M in a closed form QE.9Q . The two-dimensional integration 
over x, x' in ( |3.12| )) (and respectively, in Eqs. ( p.l| )-( |TX5| )) is to be performed numerically. 

Thus, we see that the analytical integration over the azimuthal angles is feasible and 
one end up with 2- fold integrals of the type fl3.12p instead of the 4-fold ones (|3.1|). 

For other models of MEC when the form of the matrix elements of the current operator 
is not convenient to perform the integration analytically one can utilize a procedure like 
that elaborated in [[|TJ. One can approximate a function (which may be given on a mesh) 
by a superposition of pole terms like 



F(k 2 ) * E 



C, 



{k 2 + mf) ni 



1,2,...). 



(3.18) 



We have applied this procedure in f|T| to parameterize the solution of the Faddeev equa- 
tions for the 3N bound state. Separable nucleon-nucleon potentials have been used in 
4T| to facilitate the computations. These analytical integrations for and 0' reduce the 



4-fold integrals Q to 2-fold ones of the type ( $T2 
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Another way to reduce the integrals ( |3.6| )-( |O0| ) is to write the integrand in terms of 
the spherical functions and then apply the Clebsch-Gordan decomposition 



Y hmi (p) Y hm2 (p) = 

LM 



(2/ 1 + l)(2Z 2 + l) 



1/2 



4vr(2L + 1) 
x C(kl 2 L, 000) C(kl 2 L, mi, m 2 , M) Y LM 

for products of spherical harmonics of the same argument. 
As a result, one can get integrals in the form 



(3.19) 



dp / dp' Y£, M ,(p')g (p,p',Q)y lm (P), 



(3.20) 



where the function g is a scalar. 

One can see that R^m' (pjP']Q) ~ &wm m the coordinate system where Q ~ e x . 
Therefore, instead of Eq.( |3.12| ) one can have 



-i +i 



J dx J dx' P^{x)P^{x')R M {p,p',Q;x,x'). 



(3.21) 



Certainly, this way has much in common with the one described above and in [Appendix 



D] We have chosen the first of them for implementation in the numerical calculations. 



4. Numerical Benchmarks 

We have introduced three schemes to obtain the matrix elements I 2 in (|2.18| ), e.g., the 
PWD from Appendix B, the 4FI based on ( P^B|) and the 2FI via (ggT}). In this section 
we demonstrate that each scheme works well and the results agree perfectly. We choose 
as a test case the 7r-meson exchange two-body current ( 2.24 )— ( 2.26|) . The parameters we 



use are given in Table 1. Further we choose Q=l fm _1 and put FY = 1 for the sake 
of simplicity. The PWD requires the double integrations in x and y according to ( B.2| ) 
and ( [B.4| ). Each integral is evaluated using 30 Gauss-Legendre points. The methods 4FI 
and 2FI are evaluated by using for each integral variable 10 Gauss-Legendre points. Also 
we drop the isospin part totally. Thus the quantity to be evaluated is for the spherical 
(■-component : 

J 2 (p,j/,g;(ZV)jV, (IsW) = E C(lsj, fi -m,m)C(l , s'f, f jf -m',m') 

m,m' 



X 



kn 



dpdp'Y^,_ m ,(p>) [O^p^Y^p) (s'm' \{a(2) ® a(3)}J sm) 



(4.1) 

We list in Table 2 I 2 for typical two-body channel contributions and two sets of p' and p 
values. The three schemes produce the same numbers within the 5 digits given. Note the 
change required in isospin due to ( |2.35| ). Clearly one can not deduce any physics from 
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Table 1: Parameters for the pionic current. 
he [MeV/fm] m n [MeV] l/m N [MeVfm 2 ] ff NN [MeV] 
197.33 139. 41.467 0.081 1200. 



that table, but it should be very useful for other groups for the purpose of checking their 
codes. 

To illustrate the behaviour of the currents we show in Figs. 2-6 examples for I2 for the 
seagull and pionic parts and their sum. We choose again typical channel combinations and 
threshold behaviors ~ p p can be seen. Please note the different shapes of the currents 
and the different magnitudes of the pionic and seagull parts. 

5. Summary and Discussion 

In this paper we develop three approaches (PWD, 2FI and 4FI) for calculating the matrix 
elements of the two-nucleon EM currents. We use model ||42|| , |[43|| of MEC and take into 
account the contributions originating form 7T- and p-exchanges. The considered MEC 
satisfy the continuity equation with the one-pion-one-rho-meson exchange potential. 

To get the matrix elements of interest we apply the current operator to the 3N bound 
state and project this vector on 3N basis states \pqa). The latter are specified by the 
magnitudes of the Jacobi momenta p, q and the discrete quantum numbers a in jj- 
coupling. This separation in Eqs.( [2.18| )-( |2.19| ) of both the operation and the projection 



can in general work on N-nucleon systems. 

The matrix elements in question are contained as building blocks in the amplitudes of 
the 3 He breakup by photons and electrons. We choose this particular form of the matrix 
elements since it is convenient for incorporation of the corresponding results into the 
calculations of the reaction amplitudes with full inclusion of the final state interaction. 

We express the matrix elements in terms of the multiple integrals. Resorting to the 
technique of the partial wave decompositions we reduce the matrix elements and arrive 
at the expressions convenient for the numerical calculations. 

For example, the matrix elements of MEC in \pqa) basis have the form of a superpo- 
sition of four-dimensional angular integrals. We elaborate two different ways (PWD and 
2FI) to diminish the dimensionality of integration and end up with two-fold integrals. 
One of them, 2FI, is outlined in Section 3. Expressions coming from PWD are given 
explicitly in Appendices B and C for 7T- and p-MEC, respectively. 

Performing the integration over the azimuthal angles in the closed forms (see Section 
3), we face tedious analytical calculations. To carry out the second of these integrations 
one needs to exploit the explicit form of the MEC. This means, in particular, that each 
model of MEC requires a special treatment within the 2FI procedure. 

Aiming to provide a flexibility of the approach, we work out a third procedure (4FI) 
to compute the discussed four-fold integrals with numerical methods. 

We would like to note, that the feasibility of the direct numerical calculation of the 
reaction amplitudes represented in a form of the multidimensional integrals has been 
demonstrated in [f?| , pflfl . The six-dimensional integrals evaluated there have much in 



common with the ones we deal in this work. In fact, in the both cases an essential element 
in the integrands is the mixed tensor O kK (p2,P3) arising from decomposition of the two- 
body current. The Korobov method [09 of numerical integration in many dimensions 
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Table 2: The benchmark quantity fl4.1|) for specific cases. 



Final state (/i') 


Initial state (/x) 


z/ ffm *1 


V [fill 1 

r L J - J " I " L J 


jseagull r f 31 
2 I" 11 J 


jpionic r f 3 n 




X Sn( 0) 
v w y 


0.5 


0.5 


-0 43878xl0~ 1 


0.12322xl0 _1 






0.5 


0.5 


-0 69538xl0 _1 


0.20792xl0 _1 


3 Difl) 

1 V / 


l Sn( 0) 

\ w y 


0.5 


0.5 


53929xl0 -2 

V 7 a .J '.J '7 «.7 /\ X V / 


53001xl0 -2 




3 Pnf 0) 
1 U V w y 


0.5 


0.5 


-0 13068xl0" 5 

V 7 • _L v/ V 7 V 7 v.7 / N _L V 7 


76770 xl0~ 2 

KJ a 1 KJ 1 1 v7 / \ V7 




3 Pn( 0) 

1 U\ w y 


0.5 


0.5 


-0.69148xl0 _1 


17510xl0 _1 


J- \ / 


3 Pnf 0) 

1 U\ w y 


0.5 


0.5 


-0.21208xl0 -2 


-0.15741xl0 _1 




l v y 


0.5 


0.5 


-0 43878 xlO -1 


0.12322xl0 _1 


3 Pn(0) 

* L) V w y 


J- V / 


0.5 


0.5 


-0.13068xl0 -5 


0.76770 xlO -2 


3 Pi(0) 

1 V / 


J- V y 


0.5 


0.5 


0.75258xl0 -2 


0.17771xl0 _1 


3 P1(1) 


x Pi f 0) 

-L V y 


0.5 


0.5 


25989xl0 _1 

V 7 a aW .7 ' 7 V 7 '.7 X \ _l_ V / 


-0.18852xl0 -2 


1 SnC0 , i 


3 Si(-l) 


0.5 


0.5 


-0 69538 xlO" 1 

V 7 a V7 '7 ■ ' '.J V7 / \ _1_ V 7 


0. 20792x10^ 




3 Si(-l) 


0.5 


0.5 


-0.69148xl0 _1 


17510xl0 _1 

KJ a _L 1 v7 _L W X N _1_ V7 


1 v w y 


3 Si(-l) 


0.5 


0.5 


63478xl0 _1 


-0 77629 xlO -2 

KJ all KJ kJ / N J- V7 


3 PifD 


3 Si( 0) 


0.5 


0.5 


0.54138xl0 _1 


-0.20495xl0 _1 


ip. (0) 

1 V / 


3 P1(-1) 

J- V y 


0.5 


0.5 


-0. 25989x10^ 


0.18852xl0 -2 




1 1\ w y 


0.5 


0.5 


-0 75258 xlO" 2 


-0 17771xl0 _1 

V7 a_LI 1 1 _1_ / \ _1_ KJ 




l V y 


0.5 


0.5 


-0.54138xl0 _1 


0. 20495x10^ 




3 Pi( 0) 

l v y 


0.5 


0.5 


-0 63478 xlO" 1 


77629 xlO -2 

KJ a 1 1 KJ \J / \ J_ V7 


3 Di(0) 


3 p, f-i) 

J- V y 


0.5 


0.5 


0.21553xl0 -2 


0.10227xl0 _1 


3 Difl) 


3 Pi f 0) 

J- V y 


0.5 


0.5 


-0 30655 xlO" 2 


0.11589xl0 _1 




^ 1 V x y 


0.5 


0.5 


53929xl0 -2 

V 7 a .7 '.7 ' 7 '.7 / \ -L V / 


53001 xl0~ 2 




J - / J. ^ x y 


0.5 


0.5 


-0.21208xl0 -2 


-0.15741xl0~ 1 


1 v w y 


1 V / 


0.5 


0.5 


30655xl0 -2 

V 7 a '.7 V 7 V7 7 .7 / N X V / 


-0.11589xl0 _1 




3 Di( 0) 


0.5 


0.5 


-0.21553xl0 -2 


-0.10227xl0 _1 




0) 

u v y 


1.0 


0.5 


-0 57838xl0 _1 


30361xl0 _1 

V 7 a ' J KJ v/ V7 -L / n _1_ V 7 


J- \ / 


u v y 


1.0 


0.5 


-0 34083 xlO" 1 

V 7 a ' 7 A V 7 v.7 ' 7 / \ _l_ V 7 


0.22421xl0~ 1 


J- \ / 


0) 

v w y 


1.0 


0.5 


12077xl0 _1 

KJ a _L V7 1 1 / \ _L V7 


0.11928xl0 -2 




3 Pnf 0) 

u v y 


1.0 


0.5 


-0.13462xl0 -5 


0.64118xl0 -2 


3 S 1 fl) 


U v w y 


1.0 


0.5 


-0 33735xl0 _1 


0.11544xl0 _1 


3 D,fl) 


3 Pnf 0) 

u v y 


1.0 


0.5 


-0.21350xl0~ 2 


-0.25089xl0 _1 


^nfO) 

u V / 


-L V / 


1.0 


0.5 


-0.19222xl0 _1 


16863xl0 _1 

V 7 a .-L V7 V 7 V7 '.J / \ _L V 7 


- 1 - U V w y 




1.0 


0.5 


-0.46181xl0 -5 


0.64118xl0 -2 


3 Pl(0l 

1 v w y 


l V / 


1.0 


0.5 


46799xl0 -2 

KJ a -L V7 1 kJ kJ / \ -L V7 


0.91221xl0 -2 


x 1 V x / 


x Pif 0) 

l V y 


1.0 


0.5 


0.18041xl0 _1 


-0.31948xl0 -2 


1 Snf0l 


3 Si(-l) 


1.0 


0.5 


-0.34083xl0 _1 


0.22421xl0 _1 


x U y y 


3 Si(-l) 


1.0 


0.5 


-0 95407xl0 _1 

V7 a '.7 J- KJ 1 / \ _1_ v7 


27836xl0~ 1 

KJ a 1 KJ ' J KJ / \ _l_ W 


3 Pi(0) 


3 Si(-l) 


1.0 


0.5 


0.84423X10" 1 


-0.19717X10" 1 


3 Pi(i) 


3 Si( 0) 


1.0 


0.5 


0.76826X10- 1 


-0.32340X10- 1 


^(0) 


3 P1(-1) 


1.0 


0.5 


-0.18037X10- 1 


0.31948xl0~ 2 


X P1(1) 


3 Pi( o) 


1.0 


0.5 


-0.46799xl0- 2 


-0.91221xl0- 2 


3 S!(0) 


3 P1(-1) 


1.0 


0.5 


-0.28138X10- 1 


0.13312X10- 1 


3 Si(l) 


3 Pi( o) 


1.0 


0.5 


-0.28423X10- 1 


0.11098X10- 1 
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have been applied in Refs.H-0 and | 4"I|| . This method proved to be an effective tool for 
computation of the six-fold integrals. 

The detailed analysis show that the results obtained with independent methods (PWD, 
2FI and 4FI) are in the perfect agreement. In Section 4 we present benchmark results for 
the 7T-MEC The momentum dependence of the matrix elements MEC in \pqa) basis is 
studied as well. 

The potentialities of the approaches elaborated are not restricted by the consideration 
of the simple two-body currents generated by the 7r- and p-exchange. Other, more 
complicated, models of the MEC consistent with the realistic NN potentials can be treated 
without any substantial modifications of the methods. 

So, we calculate the matrix elements of the two-body currents in a form providing the 
compatibility with approach |]l2]-|]l6fl . The results of the present article are to be used 
in forthcoming investigations of the quasi-elastic electron scattering off 3 He and proton- 
deuteron radiative capture with inclusion both the MEC and the final state interaction. 
As it is discussed above, these elements of the theory have to be treated simultaneously 
for ensuring the gauge independence of the reaction amplitudes and, respectively, for 
obtaining unambiguous predictions for the observables. 
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Appendix A. The partial wave representation of Eq. 

Eq. (|2~T9| ) reads 



W, g, Q; (A'~)/'M' V, (A^/M-//) = J dq'y^ M ,_^ q ') 6{q |g J 2 + s9H y ItU _^ + 1 ~q). 



We use the relations 



(A.1) 



7 A+2 



£ 2vr(-) fc y/khuM, q, Q)y%° k (q>, Q) 

k 



(A.2) 



with 



and 



, , , n , f 1 , „, J{q - yV 2 + Q 2 /9 + 2q'Qx/3) 

hk\{q ,q,Q) = I ^dxP k {x) 



Y lm (a + b)= J2 
i 



i h b h 



i ^|a + 6|^(2/ 1 + l)!(2/ 2 + l)!^ la '° j 



(A.3) 



(A.4) 
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Then assuming that the direction of Q is parallel to the z- axis we get employing 
standard steps 



h(q', q, Ql (A'§)I'M' - (\\)IM - = 5 rn ,^,„y \I(-) X+I '^ y/(2X + l)\ 
x Ea 1+ a 2 =a g Al (|Q) A2 ~^=^ ] Sfc H-) k h k (q', q, Q)C(\k\>, 000) 

x E/ V7C(A 2 */, 000) { f / £}{{ j/ j }c(/,J,J',0,fn-/*). 



Appendix B. The partial wave representation of Eq. 

Here we show a spherical component ( of 7 2 and drop the isospin part. This corresponds 
to Eq. ( [4.11 ). The photon momentum Q is assumed to point into the z-direction. We 
evaluated the seagull part (|2.25|) and the pionic one ( |2.26|) . 



6V3ir\Jisji's'j'(-l) l ' +s ' V 



1 1 ax 



x^(-l) ai d 1 I i i , 



oi 



1 I s 

2 2 6 

x ^(-l)^C(lja 2 ;C^ / ) 

1 «1 1 

a 3 s' a 2 
1 

-i 

2 



x E (-l) h2 (o^) Al 



x 



h 1 +h 2 =l 

^f[l — (— i) ai+/l2 +' r 



x 



Ev / Ac'(^i/i;ooo)c(/ 1 j"« 2 ;o /U V)( {j °* 



x E 4hC{rh 2 f 2 - 000)^(2/2 + 1)! 

/2 



x E (pTHpr- 





/a 


;) 




A 1} 




«3 




\ A, 


/z 2 r J 



J{2u l + l)\{2u 2 )\ 

J2^C(u 2 lz;000)C(l'zu 1 ;000)\ J U * ^ ) G h z 2 r f \ (B.l) 



where 



-i 



= /" dyP z (y)A h2 - f2 J dxP r (x) V7T (-Q 2 + A 2 + QAx) and A = \J p 2 + p' 2 - 2pp'y. 

(B.2) 
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where 



= 367rV/Wi'(-l) z+s 'V,M+C 



1 1 «! 

1 1 



2 2 S 
2 2 6 



x > e 



x 



£(-1)^(1^; C/V) 



1 Ctl 1 

x^(-l) a3 cf 3 { I s J 
«3 s' a 2 



«3 



X 



E i-V h %Q) h2 
^1+^2=1 



X 



x E ioQT 2 

x^^C(/iifflirf i; 000) 

dl 

f h ± h 2 1 1 

^Vd 2 C(/i 2 w 2 d 2 ;000) i«i t« 2 1 
d 2 [ di d 2 ct\ 

d 2 d\ ol\ 



x 



2>/SC(dilc;000) | ^ 



e c 



X 



■l) /l C(/ 1 i / a 2 ;0//V) 



x£fC(rd 2 / i; 000)C(rc/ 2 ;000) 



d 2 c e 
. h h r 

x E (p') u1 (p) U2 ^ 1 

v /(2ui + 1)!(2m 2 )! 



x 



^v / iC( M2 ^;000)C(/ / ^ 1 ;000){ ^ ^ 2 ^ 2 ) H^ wih , (B.3) 



= y dyP z (y)A 1+hM - h y d:rP r (:r) 
-l -l 

We used the abbreviation a = 2a + 1. 



^(±Q 2 + A 2 - QAx) - vMQ 2 + A 2 + QAx) 



2QAx 



(B.4) 
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Appendix C. The p Meson Exchange Currents 



The p-meson exchange current are given for instance in |E| and (43[. The current consists 
of four parts. 

3p(P2,P3) = Jp,l(pl,Pa) +Jp,Il(P2,P3) +3p,IIl(p2,P3) +1 P ,IV(P2,P3) 



JpAp^Ps) = i[r(2) x f(3)] 2 if(Q 2 ) g 

P2 ~ Ps 



3p,ii(P2,P3) = - i[r(2) x f(3)] 2 if(Q 2 ) 
x (v p (p3)a(2) x (a(3) x p 3 ) - u p (p 2 )a(3) x (a (2) x ft)) 

Jp,m(p 2 ,p 3 ) = - *[t(2)xt(3)] 2 F^(Q 2 ) ^ " ^ ((*(2) x ft) • (<?(3) x ft)) (ft-ft) 

P2 - Pi 

Jp,iv{P2,P3) = - i[t{2) x r(3)] 2 Fi (Q ) g \ 

P2 ~~ P3 

x (a(3) ■ (ft x ft)(<x(2) x ft) + a(2) • (ft x ft)(a(3) x ft)) . 
First we present their PWD for the spherical component ( : 

U p (p',p, Q; (l's')j'fi', (is)jfi) 
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where 
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where 
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1 1 «! 
1 1 



"1 



2 2 S 
V 2 2 6 



E«2 J 

"2 I 



1 1 1 

a 2 a± 1 
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x E (\Q) U2 

Ul+U 2 = l 



x E c { l 3 a Zi(w') 

ai 0:2 1 
xj^ h I s I j 



s' h «3 



x^^im/nooo) 

Ev^(lu 2 / 2 ;000) j Ma « 2 1 1 
/2 [ /1 /2 «i J 



X 



x 

912 



92 k 



#2 Z' h ) j g 2 gi a 2 
as / J 1 I h I' 



xJ2H-lY C(rf igi ;000)C(rf 2 g 2 ;000)l 9 f 2 ^ " 2 J [1 + (-1)^+-] 



1D1+U12 =gi 



x E (pO^rc-ir 



X 



v /(2 Wl + l)!(2 W2 )! 

(7l W 2 £?l 

Z Z' 2 



Ev / i(-l) z C(w; 2 /z;000)C(/ , z«; 1 ;000) 7 ^ 91 \ #*T\ (C7) 



where 



/" , s ij. /• r ,v p (\Q 2 + A 2 -QAx) -v p (\Q 2 + A 2 + QAx) 

= ~ j dyP z (y)A 1+u ^ j dxP r (x)^^ 20 Ax 

-1 -1 ^ 

(C.8) 

Secondly we present the form required for the 4-fold integration. The separation of 
spin and momentum parts like in fl2.28|) , is also displayed for the convenience of the reader. 
The first term j I,p has no spin-dependence and is left at it is. 

% Ip W,P,Q\(l'J)j'l>>',{ls)w) (C9) 
= E fi-m, m)C(l's'f, li' - m', m!) J dpdfiY^^^Yi^ip) 



xeE \ \ 7 KM^-C.C^W^®^)} 



ai v 



x[-^(p3)(-)^ +C (p 3 )c-, + ^(P2)(-) Ql ^ +C (p 2 )c-,], 

§"V 1 ?,Q;(W/*'(i#) (cio) 
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= J2 C ( ls 3i V-m, m)C(l's'j', li' - m', rn!) J dpdfY^^^Yi^ip) 

m,m' 

F2 ~~ P3 

xSEH^jl 1 ^ }w2)®^(3)} ai ,te®P3} ai ,_,, 

i? fi W,P,Q;(WW,{l8)jii) (en) 

= ^ C(Zsj, fi-m, m)C(l's'f, ll' - m', m!) J #t^^,_ m ,(^)Y^_ m (p) 

m,m' 

x vj^ z vM E(-2)v^{ 1 1 l)(-ir 

2 C(aia 2 l, k, C - «, C) {a(2) g> cr(3)} aiK ^ c ( lla 2, -« 2 - ^3 + C - «, ^2 + K3, C - «) 



1 aire 

K2K3 



XC(111,K2,K 3 ,K2 + K3) {(P2) K4 (P2)« 2 (P3) K3 + (p 3 )re 4 (Ps)k 2 (p 2 )« 3 } , 

with K4 = — — ft 3 + C — K - The momenta p 2 and ^3 are given in Eqs. ( |2.14| ). 

Appendix D. Functions F and G in the Case of it— Meson Ex- 
change Currents 



We have for the seagull contributions ( pT6D (pT 



G\m (n) - 8 m > )m 
-1 +1 



x -jdxjdx'P{r ; {x l )P{ n tx)V{m,n), 



-1 -1 



l/2\ 1 / 2 7 1 + r 

] dx dx'P[r'{x')Pr{x) (D.2) 



-1 -1 

x fp'P^(a;') V(m, n) — pP^{x)V{m + /i, n] 



and for the ones ( p.8|) -( ^T0| ) of the pionic current 



Tpl'm' _ r 

l/2\ 1 / 2+l 



x - , 

7T V3 



J J dx J dx' Pf {x^Pj^ix) (D.3) 
-1 -1 

x (p'Pf (x')W(m) - pPf (x)W{m + ll)) , 
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2 +1 +1 

—J dx J dx' P^'{x')P^{x) 



-i -i 



x ( P n P^(x')P^(x')W{m) 



(D.4) 



- pp'Pj 1 (x')fr (x)W(m + ii 2 ) 

- pp'Pt{x)P? 2 {x')W{m + /ii) 

+ p 2 P^{x')P^{x)W{m + Ml + /i 2 )) , 



x (V 3 Pf (x')Pr (^)^r (x')W(m) 

- p' 2 P Pi ll {x)Pl l2 {x , )Pi i3 {x')W{m + /ix) 

- vr (x')pr (xjpr (x'i w{m + /i 2 ) 

+ pyP 1 w (x)Pr(x)Pf(x')W / (m + /ii + /i 2 ) 

- */Vf 1 (x')Pf (x')Pr (a;) W(m + ^3) 

+ p'p 2 Pj 1 {x)P^ 2 {x')Pj i {x)W{m + /ii + Ms) 

+ pV Pf 1 (x')Pr (x)Pf (a;) W (m + H2 + Ms) 

- p' 3 Pr(x)pr(x)pr(x)^(m + ^ + ^ + ^i) 



Note, that formulae ( |D.l| )-( p3|) have been obtained in the coordinate system where 




with 



P^icosO) = V2^ e 



Y lm (6,4>). 



(D.6) 



Q || e z . 
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Appendix E. Integrals over the Azimuthal Angle 

The integrands of (|3lf )- (|3.10|) contain functions 

<V) = ^f^F nNN {ml + p 2 )- 1 and H(p 2 ,p 3 ) = (pj - p 2 2 )-\v{p 2 ) - v(p 3 )). 



In Appendix D we list the expressions for F and G (Eqs. (|3^ )-( |3.10| )) in terms of two- 
dimensional integrals over the products of the associated Legendre polynomials P^ix) 
and the integrals 

2tt 2tt 

V(m, n) = Jdt(2t — t) cos(mt)v (p n ), and W(m) = J dt(2t — t) cos(mt)H(p 2 ,p 3 ). 
o o 

(E.l) 

Functions ( [E.l| ) belong to the type of integral ( |3.14j ). To compute ( |E.1| ) analytically with 



v (p n ) and H(p 2 ,p 3 ) in question we write for the vector p n (|3.2|) 

p 2 n = a n + b cos t 

where 

1 



-Q 2 +p' 2 +p 2 + (-l) n Q(p'-p) - 2p'p cos6 cosO', 
—2p'p sine sine' 



and decompose 
v(p n 



f 2 
2n 2 m1 



a n + m 2 + b cost) 1 — (a n + A 2 + 6 cost) 1 
-(A 2 -m 2 )(a n + A 2 + bcost)~' 2 



H{P2,P3) 



NN i \-l 



2n 2 m 



— (a 2 + m 2 + 6cost) 1 + (a 3 + m 2 + boost) 1 
+ (a 2 + A 2 + b cost)" 1 - (a 3 + A 2 + b cost)" 1 
+ (A 2 - m£)((a 2 + A 2 + 6 cost)" 2 - (a 3 + A 2 + 6 cost) 



-2 



We get 
^(m, n) 



f 2 



nNN 



W(m) 



2n 2 m 2 



f 2 

2-K 2 m 2 



V(a n + ml, b, m, 1) - V(a n + A 2 , b, m, 1) 
-(A 2 - m 2 n )V(a n + A 2 n , b, m, 2)], 
NN (a 2 - as)" 1 

-V(a 2 + m 2 , 6, m, 1) + ^(a 3 + m 2 , 6, m, 1) 
+ V(a 2 + A 2 n , b, m, 1) - V{a 3 + A 2 n , b, m, 1) 



(E.2) 

(E.3) 
(E.4) 



(E.5) 



(E.6) 



(E.7) 



+ (A 2 - ml)(V(a 2 + A 2 n , b, m, 2) - V{a 3 + A 2 , b, m, 2)) (E 
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Integrals ( |3.17| ) can be calculated in closed form, e.g., with the help of tables |i7|| . For 
m > and n = 1 , 2 the result reads 



V(a,b,m,n) = 2n 2 — (\/ a 2 — b 2 — aj 
x (a 2 -b 2 )- n+ ^ (a + mV 
for m, one can use the relation V(m, n) = V(—m, n). 



I) 2 



n-l 



(E.9) 
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(a) 




Figure 2: The expression I2 from Eq. (|4.1| ), (a) the seagull term, (b) the pionic term and 
(c) the sum of both. The initial and final states are 1 Sq(/^=0) and 3 Si(// =1), respectively. 
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Figure 3: The same as in Fig. |]for 3 Pq(/U = 0) and 3 Si(/i' = 1). 
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